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An overview of theoretical results and experimental data on the thermodynamics, structure and dynamics of 
the heterophase glass-forming liquids is presented. The theoretical approach is based on the mesoscopic het¬ 
erophase fluctuations model (HPFM) developed within the framework of the bounded partition function ap¬ 
proach. The Fischer cluster phenomenon, glass transition, liquid-liquid transformations, parametric phase dia¬ 
gram, cooperative dynamics and fragility of the glass-forming liquids is considered. 
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1. Introduction 

structure of a glass-forming liquid and glass possesses a short-range and medium-range order (SRO 
and MRO) rather than a long-range order (LRO). Below the crystallization temperature, Tm, precautions 
have to be taken to avoid crystalhzation or a quasi-crystalline structure formation and to prevent the 
supercooled liquid state down to the glass transition. Therefore, a liquid can be transformed into amor¬ 
phous (glassy) solid only if cooling is fast enough to avoid crystallization. As a result, the liquid is non¬ 
equilibrium and unstable at the glass transition. For this reason a description of the glass transition can¬ 
not be based on the canonic Gibbs statistics. A palhative approach based on the bounded statistics can be 
formulated as follows. 

If the cooling time is much longer than the equilibration time of the liquid structure on scale ^ [let us 
denote this time by t (f)] and no significant structural correlation occurs on scales r > ^, one can consider 
the glass transition as a sequence of transformations of the structure states which are equilibrated just on 
the scales r < f. Statistical description of such a liquid can be developed if we exclude from the statistics 
the states with the correlation scale r > ^ and, on the other hand, ensure that the observation time, Tobs. is 
much longer than t (if) . In this case, the Gibbs partition function can be replaced by the bounded partition 
function which is used then to determine the free energy of the partially equilibrated liquid. Limitation of 
the phase space due to the exclusion of the states with correlation lengths r > f leads to an increase of the 
free energy of the equilibrium state. The standard Gibbs statistics restores with f ^ oo. The observation 
time limits from above the scale of the relaxation time t (f) and, consequently, the scale f, because t (if) 
increases ~ (the exponent 6 depends on the features of the relaxation kinetics). 

The spatial scale of the SRO, fsRO. is minimal among the possible correlation lengths in the liquid. Ac¬ 
cordingly, T (fsRo) is the shortest structure relaxation time because it is controlled by rearrangement of a 
comparatively small number of directly interacting molecules. The formation of longer correlations, with 
f » fsRO. which involves a large number of molecules in rearrangement and is driven by relatively weak 
multi-molecular forces, takes much longer than t (fsRo) time. The liquid or glass, which is equilibrated on 
the scale r ~ ^sro without considerable correlations on larger scales, is the minimally ordered amorphous 
state which can be considered using the bounded statistics method. For this reason, as the first step, the 
bounded partition function should be considered taking into account the states with equilibrated SRO. 

* Dedicated to Prof. C.A. Angell on the occasion of his 80* birthday. 
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It is experimentally established that the glass-forming liquids are heterophase (their structure con¬ 
sists of the mutually transforming fluid-like and solid-like substructures). Observations of the het¬ 
erophase structure of glass-forming liquids are numerous. Among significant observations of the last 
decades we should mention the formation of the Fischer cluster (fractal aggregate of the solid-like HPF 
in glass-forming organic licruids and polymers fillip), evolutive HPF in supercooled triphenyl phosphit 
observed in lllj], and others 

Many types of the SRO usually coexist in glasses. In Bernal’s mechanical model of the dense random 
packing of hard spheres, six types of the local order (Bernal’s holes) are statisticahy significant and nearly 
one third of them are non-crystalline flill . Similar results are obtained using computer simulations of 
liquids and glasses with different interatomic potentials iiiliil. 

A wide spectrum of relaxation times in glass-forming liquids is observed due to the variety of SRO 
types (see il^l2oj] and references cited). 

Since SRO is the molecular order formed due to microscopic forces, the correlation length ^sro is 
equal to or exceeds the range of direct molecular interactions. Therefore, to describe heterophase states, 
a mesoscopic theory is needed, in which molecular species of size r ~ ^sro specified by SRO are “ele¬ 
mentary” structural elements rather than molecules. These are not molecular potentials that determine 
the equhibrium states and relaxation dynamics of heterophase states but rather the parameters of het¬ 
erophase fluctuation interactions connected with molecular potentials. Evidently, the mesoscopic Hamil¬ 
tonian is more universal but less detahed than the microscopic Hamhtonian specified by molecular po¬ 
tentials. Parameters of the mesoscopic Hamiltonian can be considered as phenomenological coefficients 
with averaged out microscopic detahs of molecular interaction. 

These ideas are in the base of the heterophase fluctuation model (HPFM) which is con¬ 

sidered in sections [2}{^ and in appendixes|^and|B] It is further used whhe considering the issues of the 
thermodynamics of a liquid-glass transition and polymorphous transformations of glass-forming liquids 
and glasses induced by the SRO reordering and mutual ordering of heterophase fluctuations (section [7] 
and appendix[3- The cooperative relaxation dynamics of a hetei^hase liquid is considered within the 
framework of phenomenological model formulated in HPFM (section[8). Conclusive remarks are 

placed in section[9] 


2. Hetrophase fluctuations and the order parameter 

The heterophase fluctuation is an embryo of a foreign phase in the matrix phase. In many liquids, 
even in normal state (above the crystallization temperature, Tm), solid-like species are revealed by means 
of difractometry. The first observation of such heterophase fluctuations (HPF) was made by Stewart and 
Morrow 1^ . They have discovered sybotactic groups (transient molecular solid-like clusters possessing 
specific short-range order) in simple alcohols above Tm- 

The HPF are non-perturbative fluctuations in contrast to perturbative fluctuations of physical quan- 
titles near their equilibrium values in the homophase stat^. Theory of the heterophase states originates 
from Frenkel’s paper (3. Frenkel has coined the term “heterophase fluctuations” and explored the ther¬ 
modynamics of heterophase states of fluid and gas in the vicinity of the phase coexistence curve. Frenkel’s 
theory is applicable to all kinds of the coexisting phases (including the fluid and solid phases) far below 
the critical point. In this case, the amount of substance belonging to HPF is small, and thus Frenkel’s 
droplet model, with non-interacting nuclei of a foreign phase, properly describes the heterophase state. 

Frenkel’s theory fails in the case of strong HPF, when the fraction of molecules belonging to the 
“droplets” is large (for example, when it is near or exceeds the percolation threshold), and thus droplet- 
droplet interaction cannot be neglected. Besides, this theory was not generalized to include in its consid¬ 
eration the states with many SRO-types of the nucleating “droplets”. Both these restrictions of the Frenkel 
model are obviated in the HPFM. 

^Survey article (H is devoted to the physics of heterogeneous glass-forming liquids. 

^ Just a few of a huge number of papers devoted to this subject are cited. 

^Review is a good introduction to the physics of HPF. The role of non-crystalline solid embryos in vitrification of organic 
low-molecular substances {e.g., phenols) was discussed by Ubbelohde in fsill . 
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The HPFM is based on the statistics of the transient solid-like and fluid-like mesoscopic species (clus¬ 
ters) which are called 5-fluctuons and /-fluctuons, respectively. By definition, each fluctuon is specified 
by SRO. The minimal size of a fluctuon is equal to the SRO correlation length, ^sro- An arbitrary number 
of types of the s-fluctuons, m^l, can be included into consideration. 

To escape needless complications, let us assume that the fluctuons are uniform-sized with size ro 
and with the number of molecules per fluctuon equal to ko ~ r^. Thus, ^sro - ^o- This simplification is 
reasonable from the physical point of view because in the both states SRO is formed due to the action of 
the same microscopic forces, and the difference of the densities of a liquid and a solid usually amounts 
to just a few percent. The solid-like and fluid-like fractions consist of s- and /-fluctuons, respectively. 

Let us denote by N the total number of molecules of liquid and byAiy,7Vi,...,Afm the numbers of 
molecules belonging to /- and s-fluctuons, 

Nf + Ni + ... + Nm = N. (2.1) 

The total number of fluctuons is A/fiuct - NIko. 

The {m + l)-component order parameter of the heterophase liquid is determined as follows: 

Ni 

{c} = {Cf,ci,...,Cm], — (2.2) 

Cf + C\ + ... + Cni — Cf + Cs—l. (2.3) 

Evidently, c,- is the probability of the molecule belonging to i-thtype of fluctuons. Ng - Ni + N 2 +... + Nm - 
CgN is the number of molecules of the solid-like fraction. The spatial distribution of the fluctuons on scale 
r » ro can be described by the order parameter fields c/ (x) with mean values equal to c;. 

Let us regard the A:-th type s-fluctuons as statistically insignificant if Ck « m~^. The /-fluctuons be¬ 
come statistically insignificant if Cf « 1. The exclusion of the statistically insignificant components of 
the order parameter from consideration allows one to simplify the equations of HPFM. The statistically 
insignificant entities, when necessary, can be included into consideration as perturbations. 


3. The quasi-equilibrium glass transition and "ideal" glass 

Let us consider more in detail the formulated in Introduction conditions under which the glass tran¬ 
sition with equilibrated SRO takes place: 

1) The liquid cooling time or the observation time. Jobs, should be less than the time of crystallization, 

Tobs«TLRO, (3.1) 


tlro is the time of long-range ordering. 

2) The observation time is much longer than the time of short-range order equilibration, 

Tobs » Tsro ~ Ta . (3.2) 

Reordering of SRO due to localized cooperative rearrangement of the molecular structure is an elemen¬ 
tary a-relaxation event. Therefore, it is put tsro = (t^ is the a-relaxation time). 

The condition <3.11 limits the value of Xobs from above. The temperature-time-transformation diagram 
can be used to estimate Tlro and to outline the area on the {t, r)-plane in which the condition <3.1> is 
satisfied. 

The condition <3.21 restricts the value of Jobs from below. It implies that the SRO is equilibrated during 
the glass formation. Hence, the order parameter <2.31 is a function of P and T and depends on time t just 
because P and T depend on t. When this condition is satisfied, the glass transition can be considered as 
a sequence of quasi-equilibrium transformations of the SRO. 

Due to a dramatic increase of with the temperature decrease near Tg, the condition <3.21 can be 
satisfied just above Tg. Evidently, the condition <3.21 cannot be satisfied below the temperature Tp (Xobs) 
determined as the root of the equation 

Ta (T) I = Tobs • (3.3) 
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This is the temperature of kinetic glass transition because below Tp (Xobs) the SRO can be considered as 
“frozen”. Glass transition temperature Tg determined from the viscosity measurements or by means of 
calorimetry or dilatometry at the same thermal history is usually equal to Tp (Xobs) with good accuracy, 
i.e., Tp. 

In the limiting case, with Xobs ^ oo and Jobs « Xlro. when both conditions <3.11 and d3.2> are satis¬ 
fied, the quasi-equilibrium cooling of a liquid leads to the formation of hypothetical “ideal” glass (with 
equilibrated SRO and MRO but without any LRO). Hereinafter, the term “ideal glass” is used in this sense. 

It is worth to note that due to the condition (3.1) . the residual configurational entropy of the “ideal” 
glass is not equal to zero at T ^ 0 because any two parts of such a glass can be considered as non- 
correlated and statistically independent if the distance between them exceeds the largest correlation 
length which is finite by definition. 

In publications, the issues concerning the physical properties of equilibrium amorphous states below 
Tg are often debated. Between them, the hypothetical vanishing and non-analyticity of the configura¬ 
tional entropy, SconfCr), as a function of temperature, at a finite temperature Tr (the Kauzmann paradox) 
fiill . and Vogel-Fulcher-Tamman singularity of Ta (T) at a temperature Typi fsillisll are under discussion. 
In the Adam-Gibbs model the Kauzmann “entropy crisis” is included as an assumption which leads 
to the VFT relaxation time singularity at Tr. Thus, in the Adam-Gibbs model Typi = Tr. The values of Tr 
and Tvpt found from the fittings of data on thermodynamics and dynamics of many glass-forming liquids 
are close, Typr ~ Tr- Due to the above noted absence of the “entropy crisis” in the “ideal” glass, one can 
conclude that Tr and Typi should be considered as free parameters of the widely used phenomenological 
model 'he issue of proximity of Tr and Typr is considered and confirmed within the framework of 
HPFM i 



4. Mesoscopic free energy of the heterophase liquid 


The phenomenologic free energy of the heterophase liquid in terms of the introduced order parame¬ 
ter can be presented in the form of polynomial expansion in powers of {c; (x)}. 


G {P, T; {c(x)}) = Gl (P T) + Gy{PT). 


(4.1) 


In the summand GiiP T), just local interactions of the fields {c,- (x)} are included. 



(4.2) 


gL(x,Jir) = ^Ci(x)g“(P,r) + ^^Q(x)cfc(x)g°^(i^r) 

; ^ 7 ^ 


i ^ i,k 


+ TY,Ci (x)In Ci (x) -I- go (P, T). 


(4.3) 


g9 (P T) is independent of the order parameter free energy of f-th fluctuon; g9^ [PT)>0 is the fluctuonic 
pair interfacial free energy; z is the fluctuonic coordination number which is taken as independent of the 
fluctuon type. 

The summand Gv(P T] describes contribution of non-local (volumetric) interaction of s-fluctuons, 
which is taken in the following form 



r = IX - x' 


(4.4) 



(4.3) 


Wijir) is the pair correlation function of s-fluctuons, V is the volume, <l'(r) is the potential of pair inter¬ 
action of the s-fluctuons. This interaction, analogous to the attraction potential of colloid particles in a 
solvent, plays a significant role in states with diluted solid-hke species because it provides aggregation of 
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the s-fluctuons, leading to the Fischer cluster formation. It is taken as Yukawa potential with cutoff range 
Rq which is larger than but comparable with ro, 


= - — exp(-r/Ro]. 
r 


(4.6) 


Fluctuonic short-range correlation appears due to both local and volumetric interactions. The Ornstein- 
Zernike (OZ) equation (3 can be used to estimate the fluctuonic correlation length, It follows from 
OZ equation that far from a critical point, f fi is comparable with the correlation length of the direct cor¬ 
relation function, which, in turn, is comparable with the range of the fluctuonic pair interaction potential. 
With i?o ^ 2ro we have that ^fi=^2ro- 2^sro- As it is seen, the ordering of fluctuons causes extension of 
the molecular pair correlations beyond ro and the formation of the of molecular MRO. The hquid region 
of size ^fi with correlated fluctuons is referred to as correlated domain. 

The fact that the components of the order parameter At (x) are normalized probabilities, which can¬ 
not exceed 1, validates the presentation of G {P, T) in the form of the polynomial expansion in powers of 
{Ci (X)}. 

The connection of the phenomenological free energy <4.1M4.6) with the Gibbs free energy can be 
found using the approach formulated in It is shown that the free energy presented in terms 
of the order parameter plays the role of the efficient Hamiltonian in the Gibbs statistics and determines 
the most probable state of the system. The interplay between the mesoscopic free energy and the Gibbs 
statistics is considered in appendix|3 

5. The fluctuon-fluctuon interaction and the frustration parameter 

The physical meaning of the pair interaction coefficients of the neighboring fluid-like and solid-like 
fluctuons is clear. It is the fluid-solid interfacial free energy taking into account the geometry of the 
contacting fluctuons. 

The solid-like fraction can be considered as a mosaic composed of s-fluctuons with different SRO. The 
interfacial free energy of a pair of s-fluctuons depends on their mutual orientations. Evidently, coher¬ 
ent joints of the non-crystalhne s-fluctuons is hampered at any orientation. The interfacial free energy 
increase due to the geometric badness of the fit of contacting s-fluctuons is the structural frustration pa¬ 
rameter Because of its importance, let us consider the fluctuonic frustration parameter more in detail. 

A non-crystalline solid-like cluster grows due to the attachment of new molecules. Hence, the former 
surface molecules become the inner ones and the non-crystalline cluster structure becomes frustrated 
because not all newly formed coordination polyhedra are exactly similar to the initial polyhedron. A 
part of them can have the geometry similar to that of the initial coordination polyhedron but slightly de¬ 
formed. The occurrence of the coordination polyhedra of completely different geometry is also possible. 
Thus, if the initial coordination polyhedron has some symmetry, the newly formed coordination polyhe¬ 
dra have a violated or completely changed symmetry. Consequently, the binding energies of the attached 
molecules appear smaller than that of the inner molecule. 

A decrease of the binding energy per molecule is accompanied by an increase of the configurational 
entropy due to ambiguities of the geometrical changes of the new coordination polyhedra. 

As an example, let us consider the growth of a z-vertex coordination polyhedron in the case when 
the addition of a new coordination shell leads to the formation of z -1 new coordination polyhedra with 
similar but deformed initial coordination polyhedron while one of them has a different geometry. In this 
case, the energy of the inner z -t 1 molecules is 


Ez+l = £0 (z -I- 1) -t £def (z - 1) -I- £l = £o (z + 1) + £fmst 


(5.1) 


£o is the mean energy of the initial cluster, is the mean energy of deformation and £i is the energy 
of a molecule with the coordination polyhedron of different geometry. The last two terms in r.h.s. of 


^ For more information on the structural frustration see e.g., and references cited. The importance of the frustration param¬ 
eter at glass transition was considered and discussed quahtatively in . A specific frustration parameter avoiding the critical 

point is introduced in the model of frustration-limited domains (FLD) 
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O) determine the frustration energy, efrust- Because of uncertainty of the last molecule position, the 
frustration configurational entropy due to this uncertainty is as follows: 

Sfrust = Sz+I = (5.2) 

The frustration free energy is as follows: 

S^frust — £frust ~ Tifrust • (5.3) 

As it is seen, efrust is ~ z while s&ust ~ Inz. Therefore, gfrust > 0 with z » 1. 

One can conclude that generally the structure of interfacial layer of contacting fluctuons is frustrated 
and that g^ust > 0. 


6. Equations of the liquid state equilibrium 


Variation of the free energy functional (3.1) at condition (2.3) yields the equations of equilibrium state. 


S 

Sciix) 


G{P, T) + X 


d 

dci (x) 


E 


Cfc(x) = 0, 


( 6 . 1 ) 


A is the Lagrange multiplier. 

Let us denote by q; {P, T) the derivative 

(iJ T) = = g? + ^ ck {x)gik+T In a (X) 

+ I ®(x,x')Cj(x')d^x'. 
j 

Here, g/j; = Variables [P T) are not shown. 

As a result, it follows from (6.11 that 

fiflP D = Mi(f! T) = ... = ^miP T) = -A. 

These equations are analogous to the Gibbs equations of the equilibrium of phases. 
Equilibrium state is stable if the quadratic form || || is positively definite. 


( 6 . 2 ) 


(6.3) 


7. Solutions of the equations of state 

7.1. Two-State approximation 

In the physics of glass-forming liquids, different two-state models are in use for a long time 
HPFM in the two-state approximation provides abbreviated entry of the glass transition. 

In fact, in the two-state approximation of the HPFM, the mesoscopic substructure of the solid-hke 
fraction is neglected and the order parameter in the two-state approximation has just two components, 
Cs and Cf, 

Cs-tc/=l. (7.1) 

Applying the spatial averaging, we obtain from (6.2> - (6.3) 

[l-2cs]gsf+T\n-^^h,f. (7.2) 

1 Cs 

Here, 

gsf = gsf - gssl2; hsf = g° - g° - gssl2, (7.3) 

g° = Ecfcgfc+rEcpnc*, gss^Y.SikC*cl, c* = cilcs, (7.4) 

k k 
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gss is the frustration parameter. It depends on the interaction coefficients of the s-fluctuons and proba¬ 
bilities {c*\. For a while, the volumetric interactions <4.6> are not accounted for. 

In the two-state approximation, the coefficient gsf and the frustration parameter gss are taken as con¬ 
stants. Some remarks concerning the accuracy of two-state approximation of HPFM appear in section[9] 
Equation <7.21 is isomorphic to the equation of state of the Ising model with an external field hsf. The 
solution of equation <7.21 at Cs « 1 is as follows: 

CsiT) = exp{[As^, (r°) {T° -T)- gsf] /r°}. (7.5) 

Here, 

d[g°s-g°f] 

- ST ^A^T]-Ss(T] (7.6) 

is the difference of entropies of the /- and 5 -fluctuon. Tg is the solution of the equation 

8f[PA°]-gs[PA°]. (7.7) 

At Cy- = 1 - C 5 « 1 

Cf{T) = exp{ (ri) [Ag -T)- gsf] /Ag}, (7.8) 

where Tg is the solution of the equation 

gf[PTl)^gUPTl) + gss. (7.9) 

The physical meaning of the characteristic temperatures Tg, Tg is explained below. 

The temperature Tg, at which the “external field” hsf is equal to zero, is the coexistence temperature 
of two heterophase liquid states determined by equation 


g}{P,Tg]^gUPTe) + gss/2. 

At T - Tg, we have CsiTg] - Cf{Tg) - 1/2. In the vicinity of Tg, 

2Thlj,m 

^ 2 ^ 2{2Tg - gsf) ^ 3(2Tg - gsf)^ 

As it follows from <7.71 . <7.91 and <7.101 . 

Tg ~ Tg + gss/2Asf s, Tg ~ Tg — gss/2Asf s ■ 


1 ^•S 5 ,/(Te)(r Tg) ,, 

= - + — - -+ o((r- Tg)^]. 

2 2[2Tg-gsf) ’ 


(7.10) 


(7.11) 


(7.12) 


The solution <7.111 is stable at gsf [P, Tg) < 2Tg. If gsf [P, Tg) > 2Tg, it is unstable and at T = Tg, (P) the first 
order phase transition takes place. 



Figure 1. (Color online) The soUd-like fraction of of liquid, Cs, vs T/Tg at (a) gsfiPTg] > 2Tg and (b) 
gsfiP,Tg]<2Tg. 
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Figure 2. (Color online) Schematic representation of the hetrophase liquid states: (a) rare s-fluctuons 
in fluid described by equation (7.S) : (b) /-fluctuons in glass [equation (7.8H : (c) heterophase state with 
comparable fractions of the s- and /-fluctuons [equation (7.11H . 

Graphic representation of solutions of equation l l7.2t is shown in figure [1] The stable and unstable 



(7.13) 


then, the 2nd order phase transition takes place atT - Te{P). In accordance with <7.51 and <7.81 . above Tg 
and below Tg, the HPF are weak but within the temperature range [Tg, Tg], where in compliance with 
<7.111 Cs and Cf are comparable quantities, they are strong. 

It worth to note that the solutions <7.51 and <7.81 reproduce the results of Frenkel’s model in the vicin¬ 
ity of the phase coexistence temperatures (here, Tg and T^, respectively). Therefore, Tg can be consid¬ 
ered as the coexistence temperature of the fluid and heterophase liquid phases while Tg is the phase 
coexistence temperature of the “ideal” glass (as it is determined above) and the heterophase liquid. Thus, 
Tg is the ideal glass transition temperature. The real glass transition temperature, Tg, which depends on 
Tobs (see section!^, is above Tg due to dramatic retarding of the structure relaxation with temperature 
decrease. For this reason, the real glass transition temperature range, [ Tg, ], is narrower than [Tg,Tg]. 
The structure of the heterophase states in the vicinity of the characteristic temperatures Tg, Tg and Tg 
is schematically presented in figure [2] In figure [S] the mesoscopic structure of the solid-like fraction with 
several types of s-fluctuons is shown schematically. Let us remind that the solutions of equations <7.51 . 
<7.81 . <7.111 are obtained under the assumption that the fractions of s-fluctuons {c/}, / = 1,..., m are nearly 
constant or they are changing continuously and smoothly. This assumption fails if a phase transformation 
with stepwise changes of the fractions {c;} within the solid-like fraction takes place. In the next section, 
the impact of such a phase transformation within the solid-like fraction on the features of the fluid-solid 



Figure 3. (Color online) The same as in figure 2 (c) but the mesoscopic structure of the solid-like fraction 
containing several types of the s-fluctuons is shown. 
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phase transformation is considered. 

7.2. Phase transition in the solid-like fraction 

Evidently, a phase transition in the solid-hke fraction causes a non-anal 5 dic behaviour of the solutions 
of equation I l7.2h This t 5 rpe of the liquid-liquid transition appears due to multiplicity and interaction of 
the s-fluctuons which leads to the mutual ordering and phase separations within the solid-like fraction. 

As a minimal model, let us consider the heterophase liquid with two types of s-fluctuons. Hence, 
m-2. Thus, in <6.3> i,]' - 1,2. The equation of state I I6.3> for the solid-like fraction is as follows: 

c.* 

(l-2Ci)Csgl2+ Jrin—^ = hl2, C* = CilCs, Ci-i-C 2 =l, 

i Cj 

gl2 = gl2-(gl2 + g22)/2, ^12 = g2-gl + Cs(g22-gll)/2. (7.14) 

It is seen that this equation is isomorphic to equation <7.21 but the “external field” hu and the pair 
interaction coefficient Csgu depend on q. Therefore, associated solutions of equations <7.21 and <7.141 
should be considered together. The search for a general solution of these nonlinear equations at an arbi¬ 
trary set of coefficients is a cumbersome and hardly attractive task because the values of the coefficients 
for substances are initially unknown. Nevertheless, we can look for some “typical” solutions at a reason¬ 
able specification of the coefficients. 

As a useful example, let us consider solutions of equation <7.141 in the vicinity of the coexistence 
curve, hi2{P,T) - 0, assuming that Cs(g 22 “gii) is a negligible quantity. In this case, the coexistence 
temperature, T^, is determined by equation 

g?(P,ri2) = g2°(iiri2). (7.1S) 

It is assumed that T 12 is above the coexistence temperature Tg. A phase transformation of the solid¬ 
like fraction and induced liquid-liquid phase transition at ri 2 < Tg is considered in a. 

In the vicinity of T12 

hi2 = g? {P, T] - g2° {P, T) = (si - S2) {T- Tu) = Asi 2 (T- ri2), (7.16) 


Si and S 2 is the entropy of s-fluctuon of type 1 and 2, respectively. 

If Cs{Ti 2 )gi 2 > 2Ti2, then there are two stable solutions of equation <7.141 and, as a consequence, 
two stable solutions of equations <7.21 . <7.31 . Consideration of these solutions is contained in appendix[B] 
Their graphic representation is shown in figureH] The jump of the parameter Cj at T = T 12 is [see <B.71 in 
appendixjB) as follows: 

Acs (ri 2 ) ~ (si - S 2 ) exp {-gsfP) ■ (7.17) 




Figure 4. (Color online) First order phase transition within the solid-like fraction (a) induces the liquid- 
Uquid first order phase transition (b). 
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Figure 5. Continuous but stepwise evolution of cj' (T) shown in a) induces the appearance of the inflection 
point of Cs [T] at T = T 12 , b). 


If Cs{Ti 2 )gi 2 < 2Ti2, the fractions and parameter Cs change continuously at T = T 12 but the 

inflection points of the functions Cj (T) and Cs{T) appear at T = T 12 (figure [S). Hence, the expression 
dV.l?) estimates the bench height of the parameter Cs- 

The heat of the phase transitions is determined by equation dB.8t in appendixjB] 


AHu - CsiTu) (Si - S2 )A;o1 [Hf,+ Tu], (7.18) 

Hfs is the heat of the fluid-solid phase transition. It is taken into account here that the heat of 1 ^ 2 
solid-solid phase transition is equal to (si - S 2 )fcg ^ T 12 . 

It is worth to note that at a fixed value of Cg, equation (7.141 is isomorphic to Ising model with non¬ 
zero external field and with exchange integral Cggu which can be positive or negative. The external field 
controls the ratio Cj / while the sign of the exchange integral determines the type of mutual ordering of 
the 5-fluctuons. With Cggu > 0 (ferromagnetic interaction), 5-fluctuons of different types tend to separate. 
At Csgi 2 < 0, the “antiferromagnetic” order with alternating s-fluctuons of different types is preferable. 
In both cases, the fluctuonic SRO generates molecular medium-range order with the correlation length 
ifl ~ 2ro [in compliance with the general conclusion made in section |4| after equation (4.6) 1. 

7,3. The Fischer cluster 

Along with the above considered types of MRO appearing due to local fluctuonic interaction, there 
is a different type of the fluctuonic order with comparatively large (as It was observed, up to ~ 102 nm) 
correlation length, ^fc »<?//• h appears due to the aggregation of the 5-fluctuons under the effect of the 
volumetric gravitation potential (4.61 . The equilibrated aggregation of s-fluctuons possesses the fractal 
structure with fractal dimension, correlation length, equilibration time and relaxation dynamics depend¬ 
ing on the liquid features and temperature. This remarkable phenomenon, which was discovered and 
investigated in detail by Fischer et al. is known as the Fischer cluster. The Fischer cluster 

was visualized by observing a speckle pattern in ortho-terphenil @]. The speckle pattern fluctuates and 
rearranges very slowly, with characteristic time ~ 1 min at T = 293 K, while the n-relaxation time at 
this temperature is = 40 ns. With temperature increase, the speckle size and contrast decreases and 
at r > 340 K no speckle is seen. Schematically, the heterophase liquid structure with and without the 
Fischer cluster is shown in figure [6] It is worth to note that the Fischer cluster formation in heterophase 
liquid is not an exclusion but the rule if the Fischer cluster equilibration time, tfc, is shorter than the 
observation time, i.e., if Ta « Tpc « Tobs « tlro- The heterogeneous structure and slow structural re¬ 
laxation are observed not only in many Van der Waals molecular liquids but also In some metallic melts 
above Tm 

The Fischer cluster was originally identified using the results of the small-angle X-ray scattering on 
the density fluctuations. The conventional large-scale density fluctuations in a homophase liquid are 
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Figure 6. (Color online) Schematic fragments of the heterophase liquid with (a) and without (b) Fischer’s 
cluster. Just two tjqtes of statistically significant s-fluctuons are shown. The circle shows the size of the 
correlated domain. 


proportional to the isothermal compressibility kj and are independent of the wave vector q, 

\p[q)f\q^Q~KTT. (7.19) 

Here, q is the wave vector, p{q] is the amplitude of density fluctuations. The intensity of X-ray scattering 
on the density fluctuations, l{q], is proportional to |p(t 7 )p. 

It appears that a ^-dependent excess scattering intensity, /exc (^) ~ ^7”^ (D is the fractal dimension) 
occurs at T < Ta~ Tg. The /exc (^) is much larger than the scattering intensity on the thermal fluctuations 
<7.191 . The results of the wide-angle X-ray scattering show that SRO of the liquid contains both the fluid¬ 
like and solid-like components at T <Ta It turns out that the thermodynamics and a-relaxation 

dynamics are quite the same in the liquid states with and without the Fischer cluster. It means that the 
changes of the thermodynamic properties due to the Fischer cluster formation are too small to be reliably 
detected and that the fluctuonic SRO does not undergo noticeable changes. Therefore, the Fischer cluster 
formation can be considered as the process of self-organization of the correlated domains (CDs), i.e., 
entities possessing the fluctuonic SRO with the correlation length ~ 2ro, (section|4). 

Theory of the Fischer cluster is developed in (lHl23l - l^ . The Fischer cluster is considered as a fractal 
aggregation of s-fluctuons with the fractal dimension Df and correlation length ^fc- Minimization of the 
free energy as function of q, Df and ifpc allows one to determine the equilibrium values of Dp and fpc- 
It is found (see appendix[C) that 


^vc[cs,Df)~[cs)‘^f ^ ro, Dp = 3-|ln^j lnCs,cD, 


To 


Cs,CD is the concentration of s-fluctuons within CD ( 12442611 , appendix[C) 


(7.20) 


Cs.CD - 


^ ^ /(£>/,^FC) 

^d2gO(c,)/dc2-7(Dp,^Fc) 


Cs^Cs. 


(7.21) 


The quantity l{Df,^), <C.21< . is proportional to (fiQ. Equation <7.2H shows that within the CD, concentra¬ 
tion of s-fluctuons is larger than its mean value Cg’. c^.cd - Cj -i- const (/jq- 

The liquid state with the Fischer cluster is stable (while the state without the Fischer cluster is 
metastable or unstable) at 


Cs > C5,0 < 



0.16, 


l<Df^3. 


(7.22) 


Transformation of the state without Fischer’s cluster into the state with Fischer’s cluster is a weak 
first order phase transition with the transformation heat oc cpo. 

The upper bound of the Cj-range in which the Fischer cluster exists, Cs,i = {^fiDof'Csfi [see <C.26H . 
decreases, ~ (po, with an increase of the strength of the 5-fluctuons gravitation potential (pQ. When Cg 
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Figure 7. (Color online) Parametric phase diagram on the plane T* - TITg, gg^- gssl^Sg j^Tg. 

The phase coexistence lines T = Tg\ T = Tg', T = t} and the threshold of the Fischer cluster formation 
temperature, approximately determined using equations (7.St . 17.22) are shown. The critical end 
point at g*g g exists on the line T = Tg. Bold line 1 schematically presents an evolution phase curve of the 
equilibrium system. 


approaches Cjj < 1 from below, the fractal dimension Df approaches 3 and ^fc oo. It means that at 
Cs > Ci.i the solid-like fraction consists of the connected 3-dimensional solid-like clusters of size r > ^fi. 
Thus, at Cs > C 5 ,i, the topology of the heterophase hquid equilibrated on scale r>^fi changes. 

7.4. Parametric phase diagram 

In the HPFM, the structure and phase states of the heterophase liquid are described in terms of T 
and coefficients ggf, gss, hgf. It is useful to construct yjhase diagram (the parametric phase diagram) of 
the glass-forming liquid in terms of these parameter^. The parametric phase diagram of the two-state 
approximation is determined by equations <7.7t . 17.9) . d7.10t and the equation 17.22) in combination with 
17.5) . Namely, they determine the coexistence temperatures of different states in terms of the coefficients 
gsf,gss, hgf. The quantities Tg, , Tg and Tpc are coexistence temperatures of 

1) the fluid and heterophase liquid (Tg); 

2) the heterophase liquid and “ideal” glasslTg); 

3) the fluid-like and solid-like states (Tg); 

4) the heterophase liquid with and without the Fischer cluster (Tpc). 

Introducing the scaled temperature, T* - T! Tg, and the frustration parameter g*j = gss! ^Sf^s Tg, we can 
present the relations 17.8) in a dimensionless form, 

= r; = i-g;,/ 2 , Tl*^i-g:s. ( 7 . 23 ) 

The end critical point location on the fluid-solid phase coexistence curve, T* - T* [P) is located at 

8kc = 2gsf/ [ASs,f T°) - 4r; (P) > 0. (7.24) 

The first order fluid-solid phase transition on the phase coexistence curve takes place at g*j < g^^ c- 

The parametric phase diagram depicted on the plane [T*, g*s) using relations <7.231 . <7.24) and <7.22) . 
<7.5) is shown in figure [7] As an example, here is also shown one phase trajectory which becomes non¬ 
physical below the glass transition temperature T*. Within the range 0 < g*^ < g*j ^ the first order phase 
transition takes place on the phase coexistence hne T* - T*. A weak first order phase transition takes 
place on the line T* - T^^. 

^ Tentative phase diagrams of glass-forming liquid in terms of the model coefficients are introduced in O and then in . 
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7.5. Static structure factor and the order parameter restoration 

Pair correlation function of the density fluctuations, to [q, T], of the heterophase liquid with the Fis¬ 
cher cluster is ~ q~^ at ro « q~^ « At qro ~ 1, it is a superposition of the pair correlation functions 
of fluctuons, 


^ m 

a [q, T) = Cf{T]a)f{q] + ^ Ci{T)(i)i{q) = [1 - Cs{T)](Df{q) + Cs(.T]a)s{q), 
i=l 


(7.2S) 


Y m 

(Dsiq) =— T^CitDiiq), (7.26) 

CsiTi 

0f{q], 0i{q] are Fourier transfornts of the pair correlation functions of the /- and s-fluctuons, respec¬ 
tively. The cross-correlation terms ~ CfCs are omitted in t7.25t . The quantities 0f{q), COiiq] weakly de¬ 
pend on the temperature. For this reason, the equation l l7.25> can be applied to restore the order param- 
eter Cs{T) using the structure factors (Df{q), 0s(.q) measured in the liquid, fluid and glassy states 16111 . On 
the other hand, as it is shown in fiill . Cs{T) can be restored from calorimetric data using the relation 

Cs=^[Hf{T) 

^exp m]i[Hfm-Hsm] (7.27) 

which follows from equation <4.31 . Here, Hf{T),Hs{T) are enthalpies of the fluid and glass extrapolated 
in the temperature range [Tg,Tg], and HexpiTjis the experimentally measured enthalpy of the glass¬ 
forming liquid. Comparison of the results of the order parameter restoration from the structural data, 
using equation <7.25> . and from the calorimetric data, using relation 1 17.27) . gives a good chance to check 
the reliability of the HPFM. This procedure was performed using structural and calorimetric data of salol 
l^lilll. Results are presented in figure|^by scattered symbols. Solution of the equation of state in the two- 
state approximation (subsection l7.1) . in which the experimentally measured thermodynamic parameters 
and free parameter gsf are used, is presented there by a solid line. 

Let us remind that the analytic solution describes the order parameter Cs{T) of the equilibrated sys¬ 
tem. Therefore, it noticeably deviates from the experimentally determined values Cs{T) near the glass 
transition temperature, where the liquid becomes non-equilibrium. Relations )7.2S) and (7.271 . obtained 
without the assumption that the system is equilibrated, aUow us to recover the thermal history of “true” 
(in the phenomenological sense) value of Cj. 



T, K 

Figure 8. (Colo online) The solid-like fraction of salol vs T as it is found from the analysis of the calori¬ 
metric data (triangles), and from the temperature dependence of the structure factor (circles) liiLlsill . 
Arrows indicate the temperatures Tg, Tg, Tg, t}. Line presents an analytic solution of the equation of 
state in the two-state approximation with fitting parameter g^j^. 
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8. Dynamics 


8.1. a-relaxation 

Thermally activated cooperative structural rearrangements which can involve up to ~ 10^ molecules 
(la I 20 I IssHsTIl are called a-relaxation. A large amount of the molecules are involved in the rearrange¬ 
ment due to correlations. Structural rearrangement of a fluctuon also involves rearrangements of the 
neighboring fluctuons within CD of size ^ fi. Therefore, the size of cooperatively rearranging domain is 
nearly equal to ^fi. 

The activation energy of a-relaxation. 


dim a 1 

=-/ 3 =- 

dp ^ T 


( 8 . 1 ) 


depends on the order parameter. It can be presented as an expansion in powers of the order parameter 

( 3 ® 




+ pi f + p2 ^2 , 


( 8 . 2 ) 


Above Ta, the activation energy is nearly equal to E^,.. Cooperativity of the liquid d)mamics is induced by 
the s-fluctuons interaction which becomes considerable below Ta- 

Fischer and Bakai (H] have suggested that CD can be rearranged when all the molecules therein are 
in fluid-like state with correlations destroyed on the scale ^fi. This assumption leads to the following 
expression 

= n t/t.2 + (Hf -Hs] + 0 (cf). (8.3) 

(i - Jk/ 1J 

zcD ~ fi! of is the cooperativity parameter, which is the mean number of molecules within the CD; ffy, 
is the enthalpy of liquid-like and solid-like fraction per molecule. The first term is taken in the form 
proposed for random packings of spheres in (^[^. Its denominator takes into account the decrease of 
the free volume of the fluid and the numerator is equal to the activation energy above Ta- The Kauzmann 
temperature, Tr, is a fitting parameter (see comments concerning Tr in section[3 



Figure 9. The activation energy of saloi vs the reciprocai temperature Ql. 

Enthalpies Hf (T) and Hg (T) within the temperature range [ Tg, Ta] are understood as extrapolations 
of these functions measured at T> Ta and T < Tg, respectively. 

As an example of using the equation 1 18.3) Eill . the activation energy of salol was analyzed in (ill. The 
activation energy of salol vs the reciprocal temperature is shown in figure The experimental data are 
shown by circles. The curve is a result of fitting the formula (8.31 using parameters A - 967 K; Tr = 153 K, 

® There is no reason to believe that Eac (cj) is a singular function at Cj e [0,1]. 
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Te - 257 K; Ta ~ 325 K; zcd = 32; fco = 7. The input of the first term of r.h.s of I l8.3> in the activation energy 
is nearly 10% at T - Tg. 

It is noteworthy that since the main input in the activation energy makes the term proportional to 
Cs, itac iT) has the inflection point at T =; Tg. Shekel has proposed an efficient method of analysing the 
a-dynamics to check the applicability of the Vogel-Fulcher-Tamman formula and other phenomenologic 
and empiric expressions proposed for Ta(T) l7Cll472ll . He has analyzed many molecular liquids and found 
the flex points of E^ciT) identified as Tg. Coincidence of the values of Ta and Tg extracted from the 
dynamic, calorimetric and structural data (figures [8] [9) support the adequacy of HPFM. 


8.2. Ultra-slow modes and the Fischer cluster equilibration time 


Two relaxational modes are connected with the Fischer cluster. The phase transformation of a liquid 
without the Fischer cluster into the state with the cluster is controlled by nucleation and growth of a 
new phase. Since the phase transformation heat is small, the phase equilibration time is rather large 
compared with the time of elementary cooperative rearrangement Ta - The Fischer cluster equilibration 
time is determined in a. 


(pUfl) Df-l 


(8.4) 


Rearrangements of the equilibrated Fischer cluster on the scales » q > are registered as the ultra- 
slow modes iili . The relaxational rate of the ultra-slow mode is ~ q^. It is found within the framework 
of HPFM ill that 



(8.5) 


The characteristic times and tfc are proportional to and the proportionahty coefficients are rather 
large at ^fc » ^//- As it is seen, tfc » » Ta - Relations (8.4) . (8.5) are in harmony with the experimen¬ 

tal data. 


8.3. Fragility 


The fragility parameter, m, introduced by Angell in (3.01. is an important characteristic of the glass¬ 
forming hquid dynamics near Tg. It is taken as the measure of deviation of the temperature dependence 


of Ta from the Arrhenius law. There exist strong liquids, with small fragility parameter, m ~ 10, the most 
fragile liquids, with m ~ 10^, and liquids with moderate fragility in between. The fragility parameter is 
tightly connected with the structural properties and thermodynamics of a liquid, and for this reason, 
it is widely used at analysing the glass transition and classification of liquids. Angell’s definition of this 
parameter is as follows: 


d(logTa) 

^ac 

dp 

~ TglnlO 


( 8 . 6 ) 


In HPFM, the quantity Eac is determined by equation (8.3) . As it follows from (8.3) and (8.6) . 


1 



A 

[l-TdTgf 


+ zcD [HfiTg) - HATg]] ko^cs (Tg) 


(8.7) 


Hf[Tg)-Hs[Tg)~[sf-Ss)[Tg-Tg). 

The main contribution in m gives the second term within the brackets of (8.7) . It is proportional to the 
number of molecules involved in the cooperative rearrangement, zcd ~ as well as to the difference 
of the configurational entropies of the solid-like and fluid-like species. The difference of their vibrational 
entropies is comparatively small. Since these quantities can be measured regardless of Ta, equation (8.7) 
permits to check the relevance of the HPFM predictions. For example, it was found that for salol, HPFM 
gives m =: 67 . The fragility parameter of salol estimated in (3l is equal to 63. 
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As it was noted above, (subsection l8.3) . the Hf (T) and Hg [T] within the temperature range [ Tg, Ta] 
are understood as extrapolations of these functions measured at T > Ta and T < Tg, respectively. Nat¬ 
urally, the linear or quadratic extrapolation provides an acceptable result if the function is smooth and 
the higher derivatives are small. Phase transformations in the solid-like fraction lead to stepwise changes 
of Hg (r)and, consequently, to the stepwise behavior of Ta {T). In this case, extrapolations of Ta {T), de¬ 
termined by equation <8.IK from high and low temperatures into the range [Tg, Ta] cannot be properly 
fitted. Equations <7.17K <7.18K <8.1M8.3> determine the temperature dependence of Ta {T) in this case. 

In a series of experiments with some metallic glasses 17^8Oil , the fragility parameter value deter¬ 
mined using the data on Ta {T) in the vicinity of Tg and its value recovered from the extrapolated curve 
Ta [T) measured at high temperatures are completely different. As it is revealed (3, such a behavior 
of Ta {T) of Zr-based alloy Vitreloy 4 is connected with the liquid-liquid first order phase transition. In 
others melts, a transition of this t 5 q)e is assumed. 

Equations 17.171 . <7.181 . <8.11 - 18.31 provide theoretical description of this phenomenon known as the 
fragile-to-strong liquid transformation. More in detail it is considered in fsill . 


9. Concluding remarks 

The statistical basics of HPEM include substantiation of the mesoscopic efficient Hamiltonian and the 
application of the bounded statistics method whUe describing supercooled liquid states. Solutions of the 
equations of state of HPEM ascertain interplay of the thermodynamic, structural and dynamic properties 
of the glass-forming liquids. Thus, juxtaposing the theoretical predictions with experimental data (as an 
example, see figures[8][9) permits to cross-verify the adequacy of HPEM. 

The conditions of the liquid quasi-equUibrium evolution <3.11 . <3.21 determine the applicability range 
of the bounded statistics in which the states with the crystalline order on scale ^ are excluded. On 
the other hand, the amorphous states with the fluctuonic order having large correlation length are in¬ 
cluded in the statistics and the Fischer cluster is described within the framework of HPEM. Compatibility 
of the conditions <3.11 . <3.21 with the fluctuonic order equilibration on large scales should be specified. 

The hierarchy of characteristic scales of spatial correlations (starting from the local order and molec¬ 
ular size a), a < ^sro < ifi < - is connected with the hierarchy of time scales T(a) < Ta{^fi) < 

Tusv < t(^fc) [T(a) is the molecule oscillation time within the cage]. Elementary step of the nucleation 
and growth of the crystalline embryo is the cooperative rearrangement on the spatial and time scales 
^fi and Ta, respectively. Therefore, the formation of the crystalline embryos with the size larger than 
^fl takes much longer time than the SRO equilibration time. As a result, the condition <3.21 can be re¬ 
garded as satisfied when the condition <3.11 is fulfilled. Hence, the crystalline species of size ~ f// coexist 
with non-crystalline species within the solid-like fraction of hquid and in glass. As a confirmation, direct 
observations of the structure of metallic glasses by means of a high resolution field ion microscopy and 
transmission electron microscopy (see figure 6.5 in and references cited) reveal the coexis¬ 

tence of crystalhne and non-crystalline structural species with sizes of up to a few nanometers. 

The Fischer cluster equilibration (along with the ultra-slow modes) is observable only if the crystal¬ 
lization time is much longer than t(i^fc)- The crystallization heat (which is the thermod)mamic driving 
force of the crystallization) is much larger than the heat of the Fischer cluster formation. Therefore, the 
Fischer cluster can be observed just in normal liquids and in the supercooled hquids with strongly hin¬ 
dered crystallization. 

The amount of the solid-like fraction, Cg, determines the measure of the fluid-to solid transformation. 
Due to the definitive role of Cg{T), its description is an important issue of the theory. The two-state ap¬ 
proximation is a minimal model permitting to solve this problem considering the fluid-solid HPF states 
without details of the sohd-like subsystem. Evidently, this model is satisfactory if just one type of the s- 
fluctuons is statistically significant or when variations of the probabilities {c *} within the glass-transition 
temperature range are insignificant, i.e., if the mesoscopic structure of the solid-like fraction does not 
vary considerably. At the same time, estimation of the two-state approximation accuracy shows that it 
can yield acceptable results in more general cases. 

The accuracy of the two-state approximation can be estimated considering the states with transform¬ 
ing s-fluctuons. The assumption on a smooth evolution of the coefficients of equation <7.21 fails if the 
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phase transformations, similar to those considered in subsection l7.2l take place within the solid-like frac¬ 
tion. Nevertheless, even in this case, a stepwise jump of Cs{T) is comparatively small because the en¬ 
tropy jump and transformation heat at the solid-solid polymorphic transformation, as a rule, is small, 
Asi 2 ~ 10“^ while Asfg ~ 1. Therefore, one can expect that the two-state approximation is acceptable 
with accuracy to terms O ( 10 “'^ Cj) or even better. 

Figuratively speaking, in a general case, the glass and solid-like fraction of a liquid is a mosaic com¬ 
posed by the mesoscopic species of size ~ ro (figure [^. However, unless until the mutual ordering of 
s-fluctuons and the impact of the mosaic details on Cs is beyond the scope of interests, the two-state ap¬ 
proximation can be used. 

The question whether glass is a non-equilibrated highly viscous liquid or it is a non-equilibrated solid 
has rather got a conceptual sense. The thermodynamic continuity of the glass transition permits to believe 
that glass is a liquid with very high viscosity and long equilibration time. But, as a matter of fact, the glass 
near and below Tg, with Cj ^ 1, is solid with statictically insignificant amount of the fluid-like species. 
Nevertheless, it flows, like a polycrystal does, due to the diffusional-viscous flow O . Field-emission mi¬ 
croscopy of metallic glasses visualizes their grainy (polycluster) structure with sizes of grains ~ 102 nm. 
The Coble mechanism of the plastic deformation prevails near Tg in such a glass (3- The grainy 

structure of glass is the result of the existence of many centers of solidification within the liquid. There¬ 
fore, the polycluster mosaic structure of glass forms in liquids with different features of molecular forces. 
Slow structural relaxation hinders the “reclusterization” processes and the formation of “ideal” glass. 

The Fischer cluster topology changes with an increase of the sohd-like fraction. Its fractal dimension 
Df is less than 3 at q < Cs,i and it is equal to 3 at q Qn (see subsection l7.3l and appendix^. Thus, at the 
point Cs - C 5 ,i = 1 - const 00 , the topological transition takes place at which the heterophase correlated 
domains transform into homophase ones. It is important that this transition does not presuppose the 
Fischer cluster equilibration on scales fpc > ifi - This result denotes a change of the structural relaxation 
mode at glass transition considered in 18411 . 

The mesoscopic theory of thermodynamics and dynamics of the glass-forming liquids and glasses is 
connected with the microscopic approach based on the consideration of the potential energy landscape 
(see and references cited) by the landscape coarsening procedure used while deducing the efficient 
Hamiltonian (appendixl^. The coefficients of the fluctuon interaction save the memory on the micro¬ 
scopic potential energy landscape. 

A. The bounded phase space and efficient mesoscopic Hamiltonian 

Below Tm, a crystalline state is the most probable one. It occupies a phase space region Ocr of the total 
phase space Q. Q can be presented as the sum of the regions belonging to crystalline and non-crystalline 
states. 


fl - FIqJ. -j- . 


(A.1) 


Excluding Qcr. we obtain the bounded phase space belonging to non-crystalline states. The bounded par¬ 
tition function 



(A.2) 


determines the free energy of the non-crystalhne state 


G{P, T) ^-TlnZ {P,T). 


(A.3) 


In 0 ], the procedure of derivation of the equation of the free energy in terms of the order parameter 
<2.3> . G({c(x)}), (it is called the efficient Hamiltonian in (^) is expounded. It is based on the map of the 
phase space on the functional space of the order parameter {c(x)}. 



(A.4) 
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Here, a [q; x, p) are Fourier transforms of the components of the order parameter at a fixed coordinate 
(x, p) of the bounded GAT-dimensional phase space Oner- 

Performing the integration in the functional space, we have 

G(P,T)^-Tlnf exp[-G({c(^)});6]ndc;(^). (A.S) 

Selecting in this expression the contribution of long-wavelength components, q 0, and performing a 
polynomial expansion in powers of c, (x) we obtain 

GIP,T) = -rin [exp[-Heff({c});6] dci{q) + G{PT), 

qrQ«l,i 

HenPe]) = GiPT-Acix)}). (A.6) 

The term G(P T) takes into account the spatial fluctuations of the order parameter with qro ~ 1. To in¬ 
clude this summand into consideration is important in the vicinity of critical points. It generates random 
fields and has an Impact on criticality. It is shown in that at the end point on the phase coexistence 
curve, the first order phase transition can take place due to the impact of the random field. 

Within the framework of the method of cooperative variables, used while considering the gas-liquid 
critical point and systems with the Ising-type Hamiltonian 1^ . a procedure of accounting of G{P T) 
in the vicinity of the critical point is expounded. 


B. Solutions of the equations of state 


To get solutions of the equations f7.1M7.4) . 1 17.14) . 1 17.16) . let us consider the solutions of equa¬ 
tions (7.14) . (7.16) taking Cs(T] as an unknown function which should be determined later on using 
equations (7.2) - (7.4) . 

If Si > S 2 , the “external field” hi 2 is positive below ri 2 and negative above T^. The first order phase 
transition takes place and discontinuous transformation of the phase 1 into phase 2 takes place at ri 2 if 


Cs(7’l2)gl2>2ri2. 


(B.l) 


Under this condition, the solution similar to the solution (7.11) of equation (7.2) exists but it is unstable. 

To get the other two solutions near T 12 , equation (7.14) can be rewritten taking into account the 
relation (7.16) as follows: 

- 2a*Csgi2 + rin ~ 2(-Csgi2 + 2T)a* + ^T{a* f ^ Asu [T- Tu]. (B.2) 

l-2a* 3 

Here, a* = Cj a* = j - j- 

Near Tu, as it follows from (B.2) . there exist two other stable solutions. 


* 

“1,2 



Csgu-2T y'^l hi2 \ 

2T ) i Csgi2-2rj 


(B.3) 


at (a*) 2 «(c,gi 2 - 2 r) r-i. 
Equation 


C5(r)gi2 = 2r 


(B.4) 


determines critical temperatures of the system. Near a critical point, where (a*)^ s (csgi 2 -2T)T 
it follows from (B.2) 


a 


1 

2 


3A5i2(r-ri2) 

2T 


1/3 


-1 


as 


(B.5) 


Turning to the search of self-consistent solutions Cs{T), one can use the expressions found in subsec- 
tion lS.ll and (B.3) . (B.5) . At T > Te, equation (7.5) describes the required solutions Q.iCr) and Cs, 2 {T) if we 
put correspondingly 


c* = ci*i(r) = - + at(r) 


or 


c, = c 


i,2(r) = - + a*(r). 


(B.6) 
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Graphic representation of the solutions cj’ ^ (T), cj' ^ [T] and Cs,i{T], Cs, 2 {T] is shown in figurejs] 

The jump of Cj at T = T 12 is AcJ' (ri 2 ) ~ 1. Hence, in accordance with I l7.5> 

^Cs [T 12 ) - Cs,2 [Tu) - Cs,l (^ 12 ) ~ (S 1 - S 2 ) kg ^Cs (ri 2 ). (B.7) 


Here, si, S 2 are entropies of s-fluctuons. 

The heat of the phase transition is equal to 

AH = Acs{Ti2)Hfs + Cs{Ti2]Hi2-Cs{Ti2)l^-^-j^Hfs + Hi2 

= CsiTu)^^^[Hf, + Hi2], (B.8) 

ICq 

Hsf is the heat of the fluid-solid phase transition and H 12 - (^i - •S 2 ) ^ T 12 is the heat of 1 ^ 2 solid-solid 

phase transition. 


C. Thermodynamics and structure of the Fischer cluster 


The contribution of volumetric interactions into the free energy density, as seen from <4.4> . is as fol¬ 
lows: 


gv{P,T)^2n J <A>(r)Wsg(r)r^dr', r = |x-x'|, 

(C.l) 

mss(r) = (C5(x)c5(x')> = c^Wssfr) -t Cj, 

(C.2) 


where Cg is the mean value of Cs(x). 

Assuming that the s-fluctuons form fractal aggregates of dimension Df with correlation length i^pc. 
we look for the correlator u;ss(c) of the following form 


mccfr) = Cc — ( 


The condition 1 ^ Dy ^ 3 provides the topological connectivity of the Fischer cluster. It follows from l lC.2h 
<C.3t that 


(^ss — 



expC-r/^Fc). 


r>ro. 


{C.4) 


The parameters Df and fpc should be found minimizing the free energy t4.11 . 

As the first step, we find the chemical potential Ji^ {P, T) of the reference system accounting for the 
volumetric interactions of non-correlated fluctuons, at ijpc ^ 0 


T) 

(fo 


g°[P, T) + gss [P, T)Cs+T\nCs + (poCs, 
An I <I>(r)r^dr' - -AnP^cp. 


Hence, the equilibrium equation reads 


(1 


Cs 


-2cs)gsf,p+Tln- -= hsf,v 

1 Cc 


(C.S) 

{C.6) 


(C.7) 


Ssf.v = gsfiPT)- gss{P,T]/2-(pol2, 

hsf.v = g°iP,T)-g°iP,T)-gss{PT)/2-(po/2. (C.S) 

Thus, with (Uss(r) - 0, the role of volumetric interaction is reduced to a renormalization of coefficients 
of the equilibrium equation. If gsf,v < 2Te, Cg increases continuously with the temperature decrease (see 
subsection l7.1f . Just this case is considered hereinafter. 
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Let us denote by g°(Cs) the free energy of 5-fluctuon in equilibrium liquid at mss(r) = 0, and by 
gFc(cs;0/,^Fc) its value at cDssir] + 0. The difference of these quantities determines the correlation free 
energy, gcorr (cs;-D/, frc). as function of q, 

gcorr(Cs:£>/,fFc) = gTc{Cs',Df,^vc) - g°{Cs]. (C.9) 

It follows from <C.1> that 


oo 

gcorr(Cs:D/,<Fc) = 271 j <t>{r)LDss{r)r^dr' - CsAScoriT 


= -27iR^c^J 


^FC 


tCc 


^FC 


Ro + fpc 


^ Cs ^ ^’[Ro + ^vc 


Rq + ^FC 

c \2-a 
SFC 


2-a 


(p C^AScorr^ 


^0 QA5corr^» 


(C.IO) 


Ascorr is the entropy difference per 5-fluctuon due to the correlation, 


Ascorr = J({cs(r)lncs(r) +[1-C5(r)]ln[l-Cs(r)]} 

TQ 

- [C 5 InCs + (1 - Q) In (1 - Q)] ] r^dr, 

Csir] = r°r~^, = 47iCs^lc/3rQ^ = ^ {^Fc/ro]^f, 


(C.ll) 


Nf is the number of 5-fluctuons in the correlated part of the fractal; F [a, y) is an incomplete gamma- 
function, a-Df-l>0, y-rol^fi, 


ro 


3-Df 


if/J 

is the mean fraction of s-fluctuons within the correlated domain. 
Integration in IC.llt at Cj« 1 gives 


(C.12) 


As 


corr — 


3-Df 



(C.13) 


The correlation length ^fc can be estimated as follows. Noting, that in a fractal of dimension Df and 
radius ^fc 

Cs(0 = [^J , (C.14) 

we have ^ 

^Fc(cs,D/) = rote)^. (C.15) 

The fractal dimension Df is determined by relation <C.12t if equilibrium value of q.cd is known. To find 
it, let us consider the free energy of s-fluctuon within CD as a function of Cs, cd and minimize it. Denoting 
it by gcD (Cs,cd) = geo (cs: D/,^fc) . we have from O) 

gCoCCs.CD) = g°(Cs,CD) + gcorr (Cs, CD ;f3/.^Fc)- (C.16) 


Minimum of geo (Cs, cd) is attained at the value Cs, cd being the solution of the equation 


5g'CD(Cs,CD)/5Cs,CD - 0 (C.17) 

under the condition 

l3^gs,CD(Cs,CD)/5CsCD > 0- (C-18) 
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Noting that g°(Cs, cd) has a minintunt at Cj.cd = Cs and expanding tC.16> in series on degrees of Scs - 
Cs,CD - Cs, we have 


Id-'g^ics) 2 1 

gcoics.co) = - 


2 dci 


-r(2,y)( 


I [(^,+2cs6cs + 6c^,)(po 


^TC 


Rq + ^FC 


<Poc^ 


It follows from equation <C.17) - fO^ that 

Scs = 


nDf,Ocs 


where 


Thus, 


d'^g^[Cs)ldcl-I{Df,0’ 
ir(a,y) 


1 _, ^ i ^FC 


Cs,CD ■ 


1 + 


, Rq + ^FC 

nDf,^tc) 


(fo- 


(Cs)/dcf-/(D/,fFc) 


Cs- 


It is seen that 6cs ~ (po and Cs, cd S’ Cs due to the condition <C.18t . 

Equations <C.12t and l lC.22t determine the fractal dimension. Equation <C.12) gives 


D/= 3-(In y) ^ In Cs,CD, y=ro/^/;. 


{C.19) 


{C.20) 


(C.21) 


(C.22) 


(C.23) 


Since 1 ^ ^ 3, 


2s(lny) ^ In Cs,CD SO. (C.24) 

As it follows from I IC.22> and <C.24> . the fractal dimension changes within the range 1 ^ ^ 3 when Cj 

changes within the range 


Cs,0 ^ Cs ^ Cs^l, 


(C.2S) 


where 


Cs,0 


Cs,l 


1 + 




52 gO {Cs)ldc^-2I{Df,0 


1 + 


2I[Df,0 


52 gO (Cs)/5c?-27(D/,a 


-1 


> 0 , 


4cs,o < 1. 


(C.26) 


As it follows from I IC.26t . Cs,o < 0.16 with y = rg/^fi ~ 0.4, i.e., Cs,o is nearly equal to the percolation 
threshold of the solid-like fraction. The upper bound of the Cs-range in which the Eischer cluster exists, 
Cs,i, decreases ~ (po with an increase of the s-fluctuons gravitation strength cpo. It is worth to note that 
Df ^ 3 and ^fc ^ 00 when Cs approaches Cs,i from below. It means that with Cs S’ Cs,i < 1, the solid-like 
fraction consists of the connected 3-dimensional solid-like clusters of size r > ^fi. 
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reTepo(|)a3Hi craHi/i piAi/iHi/i: TepMOAi/iHaMiKa, CTpyKTypa, 

Ai/iHaMiKa 

O.C. BaKaM 

HHI 4 XapKiBCbKMM (()i3HK0-TexHi4HHM iHCTMTyT, 61108 XapKiB, VKpaiHa 

npeflcraB/ieHO or;iflflTeopeTi/i4HMx peay/ibrariB ra eKcnepuMeHia/ibHux flanux moflo TepMOflMHaMiKu, crpyKry- 
pn i flMHaMiKu rerepocfiasHnx CK/ioyTBopKDBa/ibHux piflUH. TeopeTMHHufi niflxifl Oasyeibra Ha MOfle;ii mbsocko- 
niHHHX rerepocfiasHHX cfi/iyKTyapiCi, axa 6y;ia posBi/iHyra b paMxax niflxofly o6MexeHoT craTHCTi/iHHOi c|)yHKL(M. 
Po3r/iBHyTO flBHLue K/iacrep (Piiuepa, nepexifl y c|)a3y CK/ia, nepeiBopeHHfl piflHHa-piflHHa, napaMerpwHHy c|)a- 
30By fliarpaMy, KO;ieKTHBHy flHHaMixy i cfipari/ibHiCTb CK/ioyiBoptOBa/ibHHX piflHH. 

KjiioMOBi c;ioBa: CKnoyrBoptOBanbHipiflMHM, nepexifly tpaay cmia, KJiacrep Wiinepa, no/iiMopcpi3M, 
napaMerpi/i'^Ha cpaaosa fliarpaMa 
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